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ABSTRACT 

> 

We present high resolution MERLIN radio images of multiple relativistic ejections 
. from GRS 1915+105 in 1997 October / November. The observations were made at a 

time of complex radio behaviour, corresponding to multiple optically-thin outbursts 
(N ■ and several days of rapid radio flux oscillations. This activity followed ~ 20 days 

of a plateau state of inverted-spectrum radio emission and hard, quasi-stable X-ray 
emission. The radio imaging resolved four major ejection events from the system. As 
previously reported from earlier VLA observations of the source, we observe apparent 
superluminal motions resulting from intrinsically relativistic motions of the ejecta. 
Oh' However, our measured proper motions are significantly greater than those observed 

on larger angular scales with the VLA. Under the assumption of an intrinsically sym- 
£h ' metric ejection, we can place an upper limit on the distance to GRS 1915+105 of 

^ \ 11.2 ± 0.8 kpc. Solutions for the velocities unambiguously require a higher intrinsic 

speed by about 0.1c than that derived from the earlier VLA observations, whilst the 
angle to the line-of-sight is not found to be significantly different. At a distance of 11 
kpc, we obtain solutions of v = 0.981q'q5C and 9 = 66 + 2 degrees. The jet also appears 
■ to be curved on a scale which corresponds to a period of around 7 days. 

^ ' We observe significant evolution of the linear polarisation of the approaching com- 

ponent, with large rotations in position angle and a general decrease in fractional po- 
larisation. This may be due to increasing randomisation of the magnetic field within 
the ejected component. We do not at any time detect significant linear polarisation 
from the core, including periods when the flux density from this region is dominated 
by radio oscillations. The power input into the formation of the jet is very large, 
> 10 38 erg s _1 at 11 kpc for a pair plasma. If the plasma contains a cold proton for 
each electron, then the mass outflow rate, > 10 18 g s -1 is comparable to inflow rates 
previously derived from X-ray spectral fits. 

Key words: 
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1 INTRODUCTION 

GRS 1915+105 is a distant black-hole-candidate X-ray tran- 
sient in the Galactic-plane, discovered by the WATCH in- 

* email : rpf@astro.uva.nl 



strument on board the GRANAT mission in 1992 (Castro- 
Tirado, Brandt & Lund 1992). Owing to large extinction 
(Av > 20 mag), the properties of the optical counterpart 
are unknown, but observations at radio, near-infrared and 
X-ray energies have shown GRS 1915+105 to be a highly 
unusual and energetic system. 
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Figure 1. Radio and X-ray flux monitoring spanning the sequence of radio ejections mapped by MERLIN. The upper panel shows the 
flux densities measured at 2.0 cm with the RT, and at 3.6 cm and 13.3 cm with the GBI. The middle panel shows the radio spectral index 
from simultaneous measurements at 13.3 and 3.6 cm. The lower panel shows the count rate in the 2-12 keV band by the Rossi X-ray 
Timing Explorer (RXTE) All-Sky Monitor (ASM). The optically thin state evolves towards a spectral index of —0.8 which is shown by 
the dashed line in the middle panel. 



Shortly after the discovery of the X-ray source, variable 
radio and infrared counterparts were identified (Mirabel et 
al. 1994). Subsequent mapping of the radio counterpart fol- 
lowing an outburst in 1994 revealed apparent superluminal 
motions in a two-sided ejection of synchrotron-emitting com- 
ponents from the source (Mirabel & Rodriguez 1994; here- 
after MR94). This was the first observation of superluminal 
motions in our Galaxy, the only previously measured proper 
motions of ejections from X-ray binaries (SS 433 & Cyg X- 
3) implying velocities of ~ 0.3c. Within a year, a second 
superluminal source, GRO J1655— 40, had been discovered 
(Tingay et al. 1995; Hjellming & Rupen 1995). A third pos- 
sible superluminal jet source, XTE J1748— 288, has recently 
been reported (Rupen, Hjellming & Mioduszewski 1998). 

The VLA observations of GRS 1915+105 by MR94 re- 
vealed proper motions of 17.6 ± 0.4 and 9.0 ± 0.1 mas d _1 
for two components separating from a core at position angles 
of 150 and 330° respectively. From HI absorption measure- 
ments these authors derived a most likely distance to the 
source of 12.5 ± 1.5 kpc, giving apparent transverse motions 
of 1.25c±0.15c and 0.65c±0.08c respectively. Proper motion 
and flux density ratios between the (inferred) approaching 



and receding components were consistent with simple bal- 
listic bulk motions. Solving for the proper motions at a dis- 
tance of 12.5 kpc, a true velocity of 0.92c ± 0.08c at an angle 
of 70 ± 2° to the line-of-sight was derived. Further obser- 
vations with the VLA following a radio outburst in 1995 
revealed relativistic ejection for a second time (Mirabel et 
al. 1996). However the proper motions and flux density ra- 
tios appeared to be somewhat different from those reported 
in MR94 possibly implying a change in the jet velocity or the 
angle to the line-of-sight. A summary of multiple VLA obser- 
vations of major relativistic ejections from GRS 1915+105 
is given in Rodriguez & Mirabel (1999; hereafter RM99). 

Pooley & Fender (1997; hereafter PF97) report ~ 2 yr 
of radio monitoring of GRS 1915+105 at 15 GHz with the 
Ryle Telescope (RT). They describe in detail quasi-periodic 
radio oscillations with periods typically in the range 20 - 40 
min, a phenomenon not previously seen in any other source, 
and first reported in Pooley (1995). Fender et al. (1997) 
report the discovery of infrared counterparts to the radio 
oscillations, suggesting an approximately flat synchrotron 
spectrum from cm to fxrn wavelengths. They proposed that 
each oscillation corresponds to a small ejection of material 
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Figure 2. MERLIN maps of total intensity from 10 - 12-hr ob- 
servations at the epochs illustrated on Fig. 1. The images have 
been rotated clockwise by 52 degrees to form the montage. Con- 
tour levels increase in factors of y/2 from the unit contour level 
indicated at the right-hand side of each image. Mapping fidelity is 
dominated by dynamic range considerations for the first 5 epochs, 
whilst receiver noise dominates later. 
© 0000 RAS, MNRAS 000, 000-000 



from the GRS 1915+105 system. Around the same time Bel- 
loni et al. (1997a, b) proposed that X-ray dips with similar 
periods correspond to the removal of the inner (< 200 km) 
region of the accretion disc, possibly advected into the black 
hole. The observation of a correlation between such X-ray 
dips and the rise of radio oscillations (PF97) seemed to com- 
plete a picture of repeated ejection and refill of the inner 
accretion disc, the ejecta emitting a flat synchrotron spec- 
trum from 2 cm - 2 jj,m before they cool due to adiabatic 
expansion losses. Further X-ray, radio and infrared timing 
observations (Eikenberry et al. 1998; Mirabel et al. 1998; 
Fender & Pooley 1998) appear to confirm this model, at 
least qualitatively. 

Long-term comparison of radio and 2-12 keV X-ray light 
curves in PF97 revealed an unusual period of activity in 1996 
July- August, hereafter the 'plateau' state, during which time 
the X-ray emission entered a hard, quasi-stable state, while 
the radio emission became relatively bright at 15 GHz. Ear- 
lier occurences of this state are reported in Foster et al. 
(1996) and Harmon et al. (1997). Bandyopadhyay et al. 
(1998) also discuss bright infrared emission during such pe- 
riods. It was speculated in PF97 (and earlier, in Harmon et 
al. 1997) that such plateau states may correspond to major 
radio ejections from the GRS 1915+105 system. Fig. 1 shows 
that a significant radio flare occurred around MJD 50725. 
This was the prelude to a plateau state lasting 20 days in 
which a bright, apparently optically thick, radio source was 
formed, coinciding with a quasi-stable, hard X-ray state. At 
the end of the plateau stage a second, major radio flare oc- 
curred which triggered the MERLIN mapping observations, 
the highest resolution imaging of GRS 1915+105 to date. 
Each MERLIN map epoch is indicated by a diamond in the 
upper panel of Fig. 1. 



2 OBSERVATIONS 
2.1 MERLIN 

GRS 1915+105 was observed using the Multi Element Ra- 
dio Linked Interferometry Network (MERLIN) array at 
4.994 GHz with a bandwidth of 16 MHz. MERLIN is com- 
prised of six individual antennas typically 25m in diameter 
and with a maximum separation of 217 km. Ten imaging 
runs were made as a target of opportunity program triggered 
by flux monitoring at the RT and GBI, showing the source 
had flared to around 250 mjy at 2-15 GHz. MERLIN, by 
default, samples both left hand and right hand polarisations 
allowing the full range of Stokes parameters to be imaged. 
Details of the observations can be found in the Table 1, and 
continuum and polarisation maps in Figs. 2 and 6 respec- 
tively. 

Each separate epoch includes observations of a flux and 
polarisation angle calibrator, 3C286, a point source cali- 
brator, OQ208 or 0552+398, the compact phase reference 
source, 1919+086 and the target GRS 1915+105. By nod- 
ding the array to the phase reference source, the interpolated 
phase variations due to atmospheric density variations were 
removed from the target source data. The nodding cycle of 8 
minutes on the target and 2 minutes on the reference source 
reflects the phase stability of MERLIN at this frequency. 

Phase referencing allows the registration of individual 
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Table 1. Measured quantities from MERLIN maps: angular separation (mas) from the core, flux densities (mjy), and fraction of linear 
polarisation for all well-resolved components. SE and NW correspond to the approaching and receding components of the first major 
ejection. NW1 and NW2 correspond to the apparent splitting of component NW in the later epochs (see text). CI - C3 are subsequent 
approaching ejections. M = P/I is the fractional linear polarisation, where measurable. The errors in fractional polarisation are typically 
0.01 except for the last epoch quoted where they may reach 0.04 for the weaker NW component. 



images to around 10 mas, given the changeable weather con- 
ditions during the observations. The absolute positions are 
tied to the radio reference frame (ICRF) using the calibra- 
tor 1919+086 whose position is known to around 17 mas 
relative to that frame (Patnaik et al 1992). 

Using the MERLIN d-programs, initial data editing 
and amplitude calibration were performed and the data pre- 
pared for further reduction in the NRAO's AIPS package. 
After further editing in AIPS, each dataset was run through 
the MERLIN pipeline, which images the phase reference 
source and applies the derived corrections. Instrumental po- 
larisation corrections were made using 1919+086 and the 
polarisation position angles were calibrated using 3C286, for 
which a position angle of 33 degrees for the E vector of linear 
polarisation was assumed. 

Subsequent self-calibration using the phase reference 
images as starting models was straight forward. Imaging the 
stokes parameters I, Q and U using the AIPS task IMAGR 
gives the final total intensity contour maps with superim- 
posed vectors representing linearly polarised intensity. All 
the images have been restored with a 40 mas (FWHM) cir- 
cular beam. 

All the maps are affected to some degree by the vari- 
ations in source structure and brightness during the obser- 
vations. In 12 hours, the source expands by (typically) one- 
quarter of a beamwidth, and during the first 4 epochs the 
flux-density changes were as much as 10% over each individ- 
ual epoch. Although these two effects are difficult to disen- 
tangle, and we are in the process of simulating such changes 
to quantify their magnitude, we believe that there will be no 
serious consequences for the images beyond increased uncer- 
tainties in the positions of components, at the level of a few 
mas, and the presence of sidelobes around bright compo- 
nents (at a few percent of the peak brightness) in the maps 
at early epochs. 

Relative positions of components were found by mea- 
suring radial distances of the components from a nominal 
map centre. RMS errors in this process are estimated to be 
6 mas, degrading to 15 mas for the weaker components in 
the later epochs. The error in the absolute position of each 
component is ~ 22 mas. 



A montage of all ten epochs is presented in Fig. 2. At 
least four approaching (SE, CI, C2, C3) components and 
one receding (NW) component are clearly resolved. The re- 
ceding component corresponding to ejection CI may be also 
resolved in the last three epochs (less likely, but possible, is 
a physical splitting of component NW). The maps have been 
rotated clockwise by 52 degrees, so that the mean position 
angle of the jet on the sky is 142 degrees. 



2.2 Ryle Telescope 

The source was monitored at 15 GHz using the RT at Cam- 
bridge. At least a short observation was possible on most 
days during this period. The details of the observing tech- 
nique are as given in PF97. The data, in 5-minute bins, are 
plotted in the top panel of Fig. 1. Data sampled at 32-sec 
intervals are shown in Fig. 7 for the rapid (oscillation) vari- 
ations observed at four epochs. 



2.3 Green Bank Interferometer 

GRS1915+105 was monitored at 2.25 and 8.3 GHz through- 
out 1997 using the 2-element Green Bank Interferometer. 
The observing and calibration procedures were the same as 
described by Foster et al. (1996). Random errors for the 
GBI are flux density dependent, approximately (one sigma) 
4 mjy at 2 GHz and 6 mjy at 8 GHz for fluxes < 100 mjy; 
8 mjy at 2 GHz and 25 mjy at 8 GHz for fluxes near 0.5 
Jy. We estimate that systematic errors may approach 10% 
at 2 GHz and >20% at 8 GHz. Fig. 1 displays the GBI data 
at 2 and 8 GHz during MJD 50700-50800 and exhibits both 
plateau (50730-50750) and flaring (50721-50725 and 50750- 
50756) behavior as described by Foster et al. (1996). 

It should be noted that the plateau state discussed in 
this paper is consistent with the definition in Foster et al. 
(1996) but that our definition, i.e. flat spectrum radio emis- 
sion at about 100 mjy at the same time as quasi-stable X-ray 
emission with a significant hardening of the spectrum, may 
be more specific. 



© 0000 RAS, MNRAS 000, 000-000 



MERLIN observations of GRS 1915+105 5 




750 755 760 

Date (MJD - 50000.0) 



Figure 3. Angular separation as a function of time for approach- 
ing components SE, CI, C2 and C3 and receding component NW. 
All approaching components arc consistent with unaccelerated 
proper motions of 23.6 ± 0.5 mas d -1 , significantly faster than 
those reported in MR94. The receding component, NW, has a 
proper motion of 10.0 ± 0.5 mas d _1 , and corresponds to the 
same ejection epoch as SE, on MJD 50750. Sjj' 08 . 

2.4 RXTE 

The target is monitored up to several times daily in the 2- 
12 keV band by the Rossi X-ray Timing Explorer (RXTE) 
All-Sky Monitor (ASM). See e.g. Levine et al. (1996) for 
more details. The total flux measured by individual scans is 
plotted in the bottom panels of Figs. 1 and 7. 



3 SUPERLUMINAL EJECTIONS 

Fig. 2 presents the ten MERLIN maps of GRS 1915+105. 
The images clearly show expansion of the source, with com- 
ponents on the left-hand side (south east - maps have been 
rotated clockwise by 52 degrees) appearing to move faster. 
The two components just resolved in the first image can 
be followed through all epochs as they move out and their 
brightness declines. We label these components SE and NW. 
Further components can be seen to be ejected later and are 
labelled CI, C2 and C3. The proper motions of the compo- 
nents are consistent with ballistic motions, and so we can 
extrapolate back to derive their times of formation (see Fig 
3). We estimate these to be as follows (in MJD): 

NW + SE : 50750.5l£° 8 

CI : 50753.5 ± 0.8 

C2 : 50756.6 ± 1.6 

C3 : 50758.0 ± 2.0 

The latest time of the NW+SE ejection is constrained by 
the observation of core-dominated radio oscillations by MJD 
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Figure 4. Derived intrinsic velocities and angles to the line-of- 
sight for MERLIN and VLA observations, under the assumption 
of intrinsically symmetric ejections. Solutions arc possible for the 
same angle to the line-of-sight during both sets of observations, 
but the velocity derived for the MERLIN observations is higher 
than that for the VLA observations, by A/3 ~ 0.1. 

50750.58 (Fig. 7), as well as uncertainties in the model fits. 
The individual components are unresolved in the first 4 
epochs. There is some evidence of them becoming resolved 
or breaking up into sub-components in the latest 3 epochs. 
Epochs 5, 6, 7 and 9 show that there is some emission be- 
tween the knots which may indicate the presence of an un- 
derlying continuous jet. However, given the possible prob- 
lems of varying flux density and structure through the obser- 
vations and the low signal-to-noise ratio in the later epochs, 
we cannot be certain about some of these features. The pos- 
sible bending of the jet will be discussed briefly below. 

We will concentrate our quantitative analysis on those 
components which can be followed through 3 or more 
epochs. The rapidly outflowing components are generally 
moving faster, and brighter, on the SE side. (Note that the 
apparent high brightness of the NW component in the sec- 
ond epoch is due to a blend with the core). 

From these components we can measure proper mo- 
tions, flux densities and (in some cases) polarisations. These 
we deal with quantitatively below. 

3.1 Proper motions 

An examination of the positions of the components shows 
an apparent expansion velocity at greater than 2c, thus rel- 
ativists effects are important. As a result the approaching 
components will have a higher apparent motion than that of 
the receding components. The component to the NW clearly 
has the lowest apparent speed we therefore identify it as re- 
ceding; there are 3 components to the SE which we identify 
as approaching. The data are best fit by proper motions of 

^tapp = 23.6 ± 0.5 mas d _1 
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and 



piroc = 10.0 ± 0.5 mas d 

All fits are good, with x? e d — 1> an< A illustrated in Fig. 3. 
The proper motion of 17.6 ±0.4 mas d _1 reported by MR94 
for the approaching component can be ruled out; fixing the 
proper motion to this value does not give an acceptable fit 
to the data. 

Following the method of MR94, under the assumption 
of an intrinsically symmetric ejection, we can derive 



p cos 6 = Mapp Mrcc = 0.41 ± 0.02, 



(1) 



where f3 is the velocity of the ejections expressed as a fraction 
of the speed of light, and 9 is the angle between the ejection 
and the line-of-sight. This immediately gives us a maximum 
angle to the line-of-sight (setting (3 = 1) of 9 ma i < 66 ± 2°, 
and a minimum velocity (setting cos 9 = 1) of /3 m i n > 0.41 ± 
0.02. 

A maximum distance to GRS 1915+105 can also be 
inferred for a maximum possible velocity of the ejecta of 
(3=1 



dn 



< 



ctanf? n 



Mrcc J 



/^app/^rec 

which can be expressed in convenient units as 



< 87 tan 6» 

11 



/^app Mrec 
/^app /^rcc 



kpc 



(2) 



(3) 



(for Mapp and n rcc in mas d 1 ) . The proper motions observed 
with MERLIN therefore constrain the maximum distance to 
be 11.2 ± 0.8 kpc. MR94 state that GRS 1915+105 must 
be further than the HII region G45.45+0.06 (mistyped by 
them as G45.46+0.06; Downes et al. 1980) based upon their 
HI spectral observations. While MR94 state that this HII 
region lies at ~ 8.8 kpc, in a recent study Feldt et al. (1998) 
adopt a distance of 6.6 kpc, and so GRS 1915+105 appears 
to lie at a distance of between 7-12 kpc, significantly closer 
than previously thought. Given the high observed HI column 
density to the source, from which MR94 inferred a large 
distance of 12.5 ± 1.5 kpc, we adopt a value of 11 kpc for 
the distance to GRS 1915+105. 

We can solve for the angle to line-of-sight and velocity, 
for any distance to GRS 1915+105. The angle to the line- 
of-sight 



tan 



1.16 x 10" 



/^app/-trec 
/J-app prec 



(4) 



where d is the distance to the source in kpc, and /i app and 
/i roc are in mas d _1 . Once 9 is calculated, we can immedi- 
ately derive (3 as the value of (3 cos 9 is already known. 

Table 2 lists the apparent velocities and solutions for 
the angle to the line-of-sight and intrinsic velocities for both 
these MERLIN observations and the VLA observations of 
MR94, for assumed distances of 9, 10, 11 & 12 kpc. Fig. 4 
shows these solutions for all distances between 7-12 kpc. 
Our MERLIN data unequivocally imply a higher intrinsic 
velocity, by A/3 ~ 0.1, in comparison to the observations 
reported in MR94. 




Angular separation from core (milliarcseconds) 

Figure 5. Flux densities of individually-resolved components as 
a function of angular separation from the core. Simple power-law 
and exponential decay models do not fit the data well. The flux 
density ratio of SE : NW, approaching and receding components 
from the same ejection, ranges from 6 - 10. 



3.2 Flux densities and ratios 

The ratio of flux densities from the approaching and reced- 
ing components is another important diagnostic of the ejec- 
tions. In the above calculations we have assumed geometric 
symmetry; here we also assume that the jet components on 
both sides have the same intrinsic luminosity (which does 
not appear to be the case for GRO J1655-40; Hjellming & 
Rupen 1995). When measured at equal angular separations 
from the core (i.e. the same time since ejection in the rest 
frame of the ejecta) and coupled with the observed spectral 
index we can test the theoretical predictions for Doppler 
(de-)boosting of the components. 

The mean spectral index of the ejecta is hard to de- 
termine due to both overlapping ejection events and the 
presence of rapid flat-spectrum oscillations from the core 
(PF97 and below). Our best estimate for the index, defined 
as a = A log S„/A log v is around -0.8, the same as reported 
by MR94. This is determined primarily from the simultane- 
ous Green Bank monitoring at 13.3 and 3.6 cm and should 
be a good estimate for the 6 cm MERLIN observations. It 
is difficult to measure accurately the ratio of flux densities 
of the components at equal angular separations, given their 
very different proper motions, but it lies between 6 and 10. 
This is compatible with the 8 + 1 reported in MR94. 

For bulk motions at velocity (3 the observed ratio is 
predicted to be 



5* a 



Sri 



1 + (3 cos 9 
1 - p cos 9 



(5) 



where, theoretically, k = 2 for continuous jets and 3 for 
discrete components. As /3cos# is already calculated as 
0.41 ± 0.02, and independent of distance, we can solve for 
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Table 2. Apparent velocities, solutions for the angle to the line-of-site, intrinsic velocity of ejection (assuming symmetry), and bulk 
motion Lorentz factors for values of proper motion measured by MERLIN and previously by the VLA (MR94), for assumed distances of 
9, 10, 11 and 12 kpc. The functions are plotted in full, with consideration of measurement errors, in Fig. 4 



k. For a flux density ratio of in the range 6 - 10, k = 1.3 
- 1.9. As in MR94, the flux ratio appears to be closer to 
that expected for a continuous jet than for discrete compo- 
nents (although the value derived from the data of MR94 
is k = 2.3, interpreted by them as implying something in- 
termediate between a continuous jet and discrete ejections) . 
For k = 3 we would have expected a flux ratio of ~ 27. 

Fig. 5 shows the flux density of each well-observed com- 
ponent as a function of angular separation from the core. As 
is obvious from Table 1 and the maps of Fig. 2, the ejections 
in the sequence NW+SE, CI, C2, C3 are of steadily decreas- 
ing flux density. While uncertainties in both flux density 
and position make interpretation difficult, the data do not 
appear to show the straight line behaviour expected for a 
power-law decay (for the case of adiabatic expansion losses 
and a constant expansion rate). Exponential fits to the data 
provide no improvement. RM99 discuss a steepening of the 
decay of the radio flux density at an angular separation of 
~ 1 arcsec; as we only image on smaller scales we cannot test 
this. Atoyan & Aharonian (1997) discuss the implications of 
the flux ratio in considerable detail. 

3.3 Polarisation 

Polarisation images have been made for the first four epochs, 
when the source was sufficiently bright to detect linear po- 
larisation at the level of a few percent. These are shown in 
Fig. 6, where the vectors represent the fractional linear po- 
larisation and observed electric field angle. The fractional 
linear polarisations are also listed in Table 1. No circular 
polarisation was detected at any time, with a conservative 
upper limit of 2%. 

3.3.1 Resolved ejecta 

The images of the ejected components show a striking asym- 
metry in the linear polarisation - only the approaching SE 
component appears significantly polarized. Its fractional po- 
larisation decreases rapidly from 14% to 6%, and the polar- 
isation position angle swings by approximately 75 degrees 
between the second and third epochs, then swings back by 
45 degrees. In the fourth epoch image, there is marginal de- 
tection of polarized emission in the receding NW component, 
at a level of 13 ± 5%. 

The changes in polarisation seen in the MERLIN im- 
ages could be a result of changing Faraday effects (internal 
or external) or changing magnetic field geometry within the 
radio components. With only a single frequency it is hard 



to distinguish between these possibilities. If the observed 
change in position angle is due to Faraday rotation, the im- 
plied change in rotation measure is > 300 rad m~ 2 . 

The radio components are believed to be considerably 
denser than their surroundings (observed limits on the de- 
celeration imply a density contrast of at least 10) and have 
cquipartition field strengths of order 50 mG (see below), 
presumably much greater than in the surrounding medium. 
The Faraday depth, or rotation measure, is proportional to 
nB(Ld) 1/2 , where n is the thermal electron density, B is the 
magnetic field strength, L is the path length and d is the 
field tangling scale. The Faraday depth within the source is 
therefore likely to exceed the Faraday depth of a surround- 
ing medium distributed on a scale size comparable to the 
source itself. As the components expand their internal Fara- 
day depths will fall, and any internal Faraday depolarisation 
would decrease. This is contrary to the observed decrease 
in fractional polarisation. Therefore, we do not believe the 
changes in polarisation are due to changing Faraday depths, 
and this is consistent with the low rotation measure observed 
by Rodriguez et al. (1995). 

Rodriguez et al. (1995) noted a fractional polarisation 
of approximately 2% at 5.0, 8.4 and 15 GHz at a position an- 
gle of 178 ± 10 degrees on 24 March 1994, when the source 
was barely resolved by the VLA and 5 days after the ra- 
dio components were ejected. At this time the emission was 
probably dominated by the approaching component. The 
similarity of the fractional polarisation and position angle 
at the three wavelengths implies that Faraday effects were 
small at 5.0 GHz. The implied upper limit to the Faraday 
rotation measure is 50 rad m -2 . 

Superluminal jets in AGN often show complex radio 
polarisation behaviour: rapid variations in polarisation posi- 
tion angle (e.g. 0917+624 for an extreme case, Quirrenbach 
et al. 1989) and small-scale variations along the jet (e.g. 
3C454.3, Kemball, Diamond & Pauliny-Toth 1996). This is 
usually interpreted as being due to shocks, which compress 
and order the magnetic field (Cawthorne et al. 1993), and 
which appear as radio knots in an otherwise continuous jets. 

The radio components of GRS1915+105, on the other 
hand, are thought to be dicrete plasmons, rather than fea- 
tures (shocks or otherwise) in a continuous jet. In this case 
the changes in radio polarisation indicate changes in the in- 
ternal magnetic field geometry, such as randomization of the 
magnetic field as the plasmons evolve. The effects of aber- 
ration on a planar structure such as a shock would result in 
the two components having different observed polarisation 
characteristics. However the low signal to noise for the mea- 
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Figure 6. Linear polarisation E vectors superimposed upon total intensity contour maps of the first four epochs of MERLIN observations. 
For polarisation vectors, 100 mas = 3.33 mjy/beam. This clearly shows the dcpolarisation of the core and receding components with 
respect to the approaching component, and the rotation of polarisation vectors in the approaching component. 



surement of polarisation in the receding component does not 
allow us to make such a distinction in this case. 



3.3.2 Core 

The core, sometimes blended in our images with newly 
formed major ejections, does not at any time show signifi- 
cant linear polarisation. In particular, during the first two 
epochs of observation, the core is dominated by the flat spec- 
trum oscillations of the type discussed in detail in PF97 (see 
Fig. 7). The lack of observed polarisation may be due to 
superposition of multiple components with different polar- 
isation position angles within the beam, or large Faraday 
depths close to the binary system (either internal to the 
compact base of the jet or in the form of circumstellar ma- 
terial) . 



3.4 An underlying curved jet ? 

In several of the MERLIN maps there is a hint of extended 
emission between the components. Furthermore, the low 
value of k derived from the flux density ratios and f) cos 6 is 
more consistent with a continuous jet than physically dis- 
crete components (see above). In addition, the jet appears to 
be curved, an effect which may also be marginally apparent 
in reinspection of the VLA images of MR94. 

Firstly, comparing the flux in discrete components in 
the MERLIN maps with respect to radio monitoring with 
GBI and RT (in which the source and jets will be unre- 
solved), we find that the values agree to within a few per 
cent. So the overwhelming majority of the synchrotron emis- 
sion, at cm wavelengths at least, arises in unresolved (with 
MERLIN) bright spots. On the other hand, we both derive 
(section 3.2) and apply (section 4.1) k ~ 2, the value the- 
oretically expected for a continuous jet. We must conclude 
that the data are not good enough at present to definitively 
establish whether the ejections are in discrete components 
or we are observing the bright parts of a continuous jet. In 
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the latter case we can at least now state that any underlying 
continuous jet is very faint and less than a few mjy at cm 
wavelengths 

The possible curvature of the jet is not addressed quan- 
titatively here, but we do believe this effect to be real. Dis- 
crimination between ballistic and helical effects, the latter 
possibly due to Kelvin-Helmholtz instabilities, is not possi- 
ble at this stage. The 'period' of any bending observed at 
angular scales imaged with MERLIN is around 7 days. It 
is interesting to note that the 'period' of the jet bending in 
GRO JI655-40 reported by Hjellming & Rupen (1995) of 
3.0 ± 0.2 d turned out to be very close to the subsequently 
discovered orbital period of 2.6 d (Bailyn et al. 1995). We 
are confident that the apparent bending is not due to the 
rapid angular and flux density evolution of the source duing 
a 12-hr observation (c.f. however, VLBI observations of GRO 
J1655-40, Tingay et al. 1995, where both higher proper mo- 
tions and finer angular resolution exacerbated the effect). 

3.5 Collimation and expansion of components 

From the lack of clearly resolved structure perpendicular 
to the jet axis at the 40 mas resolution of MERLIN, we can 
place constraints on the opening angle and lateral expansion 
of the ejected components. The maximum distance from the 
core that we track ejected components is ~ 300 mas, and 
this constrains the opening angle of the jet to < 8 degrees 
over distances from the core < 4000 A.U. (6 x 10 16 cm). 
Similarly the ratio of lateral expansion to forward velocities 
is constrained to be < 0.14, i.e. a maximum lateral expansion 
velocity < 0.14c (for jet bulk velocity of c). 



4 DISCUSSION 

4.1 Energetics and mass flow 

4-1.1 Estimation of the internal energies of the jet 
components 

Initially, we must work in the rest-frame of the emis- 
sion region. Converting the observed parameters (emission 
frequency, luminosity etc) requires the value of the bulk 
Lorentz factor 7, which is somewhat uncertain as discussed 
in section 3.1, on account of the uncertainty in the distance. 
Because of the large angle to the line-of-sight, the Doppler 
factors for both the approaching and receding components, 

<5a PP ,rcc = 7 _1 (1 =F/3cos6») _1 (6) 

are both less than unity for the range of Lorentz factors 
considered here. 

The estimation of the parameters of the synchrotron 
emission region is detailed by Longair (1994) and Hughes 
(1991). We adopt the formulae in the summary by Lon- 
gair (allowing for the opposite convention for the sign of a). 
These include a number of simplifying assumptions which 
are unlikely to introduce larger uncertainties than those im- 
posed by our imperfect knowledge of the source parameters. 

In order to estimate the internal energy of the emitting 
component, we need to know the synchrotron spectrum and 
the source geometry. A radio spectrum of the form S tx v a 
arises from an electron population with an energy spectrum 
N(E) tx £T (1 - 2a) . The observed luminosity can arise from 



various combinations of magnetic flux density B and particle 
density, both of which we assume to be spatially uniform; the 
total energy is minimised to derive the "minimum energy" 
conditions: the total energy in relativistic electrons, and the 
misleadingly-named -B m in, the field at which the total energy 
is minimised. As is well known, the minimum energy result 
nearly coincides with equal energies in magnetic field and 
particles, or equipartition. There are physical plausibility 
arguments for such a situation, but very few measurements 
which confirm that this situation might actually pertain in 
any particular case (see Harris, Carilli & Perley 1994 for 
data on the radio galaxy Cygnus A which suggest, in that 
case, that conditions may be near to those suggested by this 
method) . 

The calculation can easily be adapted to include a con- 
tribution to the energy by a hypothetical population of rel- 
ativistic protons, which do not radiate significantly by the 
synchrotron mechanism if their energies are comparable with 
those of the electrons. 

If the energy spectrum of the electrons is relatively steep 
(a steeper than —0.5), as in this case, the total energy of the 
particles is dominated by the lower end of the distribution. 
The lowest frequency detected, and the luminosity at that 
frequency, are then parameters in the calculations. 

The source geometry is, so far, not well-defined. Di- 
rect observations only give an upper limit to most dimen- 
sions: the resolution of MERLIN is 50 mas, corresponding to 
8 x 10 13 m, or 3 light days, at 11 kpc. We have to rely on the 
time-scale of the variations in flux density for a better limit; 
the rise times are the shorter and therefore more restrictive. 
The rise time may represent the interval during which rela- 
tivistic material is injected, or alternatively a transition in 
optical depth (or both). We believe that the former (injec- 
tion timescale) is more probably the dominant effect, since 
the spectral data suggest that the source is optically thin as 
soon as the jet is visible. 

Since the rise-time for the jet components appears to be 
less than 12h, we adopt 12 light-h (1.3 x 10 13 m) as a "typi- 
cal" size on formation of the jet component. In comparison, 
RM99 discuss a geometric mean angular size of 35 mas, cor- 
responding to a linear scale of 6 x 10 13 m. We stress the weak 
links in this argument, elaborated by Longair (1994) and by 
Hughes (1991): 

• We have only rather weak constraints on the source size, 
and certainly do not know the filling factor of the source, and 
this is one of the major uncertainties in the calculation. A 
smaller size or filling factor reduces the total energy required 
as (volume) 3 ^ 7 . 

• There is little evidence that the minimum-energy con- 
dition is achieved. 

• The contribution of protons (and other nuclei) is not 
known. 

• The limits on the energy spectrum are not well-defined 
by current observations. 

We now apply the formulae from Longair (1994); his 
equations 19.29 and 19.30, using the following measured, 
estimated or derived parameters: 

(i) adopted distance 11 kpc 

(ii) 7 = 5.0 (see Table 2) 

(iii) /3cos<9 = 0.41 
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Table 3. A summary of the derived properties of the ejecta at minimum energy conditions, calculated for a distance of 11 kpc, and 
summing both sides of the ejection. We assume that radio emission observed by GBI and RT at the peak of the flare are dominated by 
the approaching (SE) component, and take 1 and 15 GHz as upper and lower bounds of the observed flaring emission respectively. We 
assume a generation time for the ejection of 12 hr. Proton mass and mass flow rate are for the case of 1 proton for each electron. 



Minimum energy 


condition 




Ejecta 




Power 


Mass outlflow 


E m in (erg) E Kl min 


B min (mG) 


N e 


M e (rcst) (g) 


M P (g) 


(erg s" 1 ) 


rate (g s _1 ) 


2 x 10 43 2 x 10 44 


280 


3 x 10 46 


3 x 10 19 


5 x 10 22 


2 x 10 39 


10 18 



(iv) spectral index a = —0.8 

(v) Doppler factors c5 ap p,rcc = 7 _1 (1 T ficosOy 1 = 
0.34,0.14 

(vi) k — 2 (section 3.2) 

(vii) ratio of relativistic proton energy to electron energy, 
(77 - 1) = 

(viii) total source volume V = 10 39 m 3 

(ix) monochromatic luminosity L v in the rest frame, at 
the lowest frequency detected, as derived below. 

The maximum observed flux density observed at the 
start of the outburst at 2.3 GHz is 550 mjy (Fig 7), some 
2 days before the first MERLIN map. This would corre- 
spond to Uy at 1 GHz, assuming the spectral index to be 
—0.8. We adopt 1 GHz as the lowest frequency at which 
emission was detected during this outburst (Hannikainen & 
Hunstead, private communication, detected the outburst at 
843 MHz). We assume that this flux density is dominated 
by the approaching component. Its apparent flux density if 
observed in its rest frame would be about 20 Jy; signals re- 
ceived at 1 GHz would have been emitted at v — 3 GHz, and 
the luminosity L v = 2.9 x 10 17 WHz -1 . 

The minimum total energy is 

W min ~ 3.0xl0 6 i ) 4/7 (V/m 3 ) 3/ >/Hz) 2/7 (L„/WHz-') 4/7 J(7) 
and the associated magnetic field is 

B min ~ 1.8( J? (L„/WHz- 1 )/(l//m 3 )) 2/7 (;,/Hz) 1/7 T (8) 

This set of parameters, for one of the radio-emitting 
components, leads to VK min = 7.7 x 10 35 J = 7.7 x 10 42 erg 
and B min = 28 ^T = 280 mG. 

The Lorentz factor of an electron radiating near 3 GHz 
is then about 90, and the mean Lorentz factor for this pop- 
ulation is 240. (Some authors consider electron distributions 
which continue to much lower energies, but as yet there is no 
observational evidence for this. RM99, on the other hand, 
adopt a higher mean Lorentz factor of 1000, on the basis 
of a detection of the event at 240 GHz. That emission may 
well have originated in the inner region of the system; Fig 
7 shows the quasi-periodic oscillations observed at 15 GHz 
during the early stage of the current outburst which have a 
flat spectrum (Fender & Pooley 1998) and, with a time-scale 
of 20min, must come from a very small region.) 

4-1.2 Total number of relativistic electrons 

By integrating the electron distribution we can derive an 
estimate of the total number of relativistic electrons (this 
may become important if we suppose that each is accompa- 
nied by a "cold" proton) . This integration requires minimum 
and maximum energies, but as in the derivation of the total 
energy the precise value of the upper limit is unimportant; 



the number is dominated by the low-energy electrons. To 
establish the constants involved we use equation 19.17 from 
Longair, which relates the luminosity, the magnetic field and 
the spectrum, and derive 

iVtotai = (L,/WHz- 1 )/(a)/(B min /T) 

where /(a) = 1.7 x 10 24 for a = -0.8. This results in 
iVtotal = 1.6 x 10 46 . 

4-1.3 Kinetic energy 

Associated with the bulk motion of the jet there is also the 
kinetic energy; if the jet contains only e + e~ plasma and 
magnetic field, the magnitude of this is (7 — 1) x W m i n , 
perhaps 9 x 10 36 J = 9 x 10 43 erg. On the other hand, if 
the plasma contains cold protons, the kinetic energy will 
be increased by (7 — 1) x iVtotai ?%> c 2 . This results in an 
increase of only a factor 2. The total of the kinetic energy 
of the protons could vary, either way, if the conditions were 
far from "minimum energy", or be substantially larger if 
the electron energy distribution continues to rise to lower 
energies, with correspondingly more protons. 

We can estimate the power required to produce this in- 
jection of energy over 12 h; including the kinetic energy of 
the system and both sides of the source, we have to gen- 
erate 2 x jWmin in that time, a power of 2 x 10 32 W or 
2 x 10 39 ergs _1 . Including a similar number of protons ap- 
proximately doubles that total and also requires a mass- 
flow rate of > 10 18 gs _1 . X-ray spectral fits (e.g. Belloni 
et al. 1997b) have been used to derive accretion rates near 
10 18 gs _1 . 

Many of the observed parameters used here are uncer- 
tain, with the dominant effect arising from the unknown 
distance. Were we to adopt a distance of 9 kpc, the bulk 
Lorentz factor would be 1.8 rather than 5 and the Doppler 
factor of the approaching component near 1. Therefore its 
flux density in the rest frame would be 1 Jy rather than 
20 Jy, and the energy of the system reduced by a factor of 
about 10. 

The results of these calculations of energetics and mass 
flow are summarised in table 3 (for a distance of 11 kpc). 
Indeed the mass of cold protons could be as much again 
as in table 3 if we use the upper limit on electron density 
(< 2 cm~ 3 ) inferred from Faraday effects discussed in section 
3.3.1. 

4.2 Comparison with previous major ejections 

While it could already be guessed from the higher observed 
proper motions, our solutions (Table 2, Fig. 4, under the 
assumption of intrinsic symmetry) illustrate that we cannot 
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avoid deriving a higher velocity for the ejecta than derived 
by MR94. Thus the two data sets are not compatible with 
a simple change of the angle to the line-of-sight of a jet 
which is intrinsically physically identical. In fact the data 
are marginally consistent with the same angle to the line 
of sight for the jet. The increase in intrinsic (3 from our 
observations compared to MR94 is at least 0.05, and more 
likely ~ 0.1, corresponding to a significantly higher Lorentz 
factor for bulk motion by a factor of > 1.3, as tabulated in 
Table 2. 

It is of great importance to the determine whether the 
different measured velocities are intrinsic to the ejecta or 
an artefact of the different resolutions of MERLIN and the 
VLA. In the former case the significantly different measured 
proper motions could correspond to a difference between 
the ejection event (s) reported here and the event reported 
in MR94. The earlier event appears to have been both sig- 
nificantly brighter and to have decayed in flux more slowly. 
Alternatively, there may be genuine deceleration between 
the angular scales of up to 300 mas imaged with MERLIN 
and greater than 400 mas measured with the VLA. This 
does not seem likely however, given that both instruments 
recorded essentially ballistic motions (ignoring the slight ap- 
parent curvature) within their multiple observations. 

It may be that the different measured proper motions 
arise from an inability of the VLA to resolve individual com- 
ponents, causing them to blend together in maps. In this 
case, the VLA could measure lower proper motions and/or 
apparent decelerations if 

(1) There were multiple components, and 

(2) The components decreased in flux density more 
rapidly with distance from the core. 

This effect is discussed by Hjellming & Rupen (1995) 
for GRO J 1655-40 where they measure apparent deceler- 
ations of ~ 30% with the VLA in comparison to VLBA 
(and SHEVE - Tingay et al. 1995) for two of the three ma- 
jor ejections that they image. Did these effects occur during 
the observations of MR94 ? Their maps do indeed show that 
multiple ejections were occurring during this period (MR94; 
RM99); and furthermore RM99 discuss evidence that the 
decay rate of the ejecta increase with angular separation 
from the core (although they discuss a relatively abrupt in- 
crease in decay rate at angular separations > 1 arcsec, which 
would not have affected the measurements of MR94) . So it 
seems a possibility that the lower measured proper motions 
arise from blending of multiple components and that the 
MERLIN measurements more accurately represent the true 
situation. 

In order to test this we have attempted to convolve our 
MERLIN data with lower-resolution beams comparable to 
that of the VLA observations of MR94. Note that the com- 
parison is not good because of the lack of short baselines 
in the MERLIN array compared to the VLA. Nevertheless, 
it was clear from convolution of our data with a 250 mas 
beam that it was not possible to clearly resolve any individ- 
ual components. Even with a 100 mas beam, i.e. twice the 
resolution of the VLA data of MR94, while the approaching 
and receding sides can be distinguished, individual ejecta on 
the approaching side cannot. 

However, without simultaneous VLA and MERLIN 
(and possibly also VLBI) observations of the outburst, it 



is not possible to test this conclusively. So we cannot con- 
fidently discriminate between the two most likely expla- 
nations, of either intrinsically different velocities between 
March/April 1994 and Oct/Nov 1997, or resolution effects 
between the two arrays. Note that if the latter is the cause 
of the different measured proper motions then it requires the 
similar proper motions (range 15 - 18 mas d _1 ) measured 
in four instances in 1994 with the VLA (RM99) all to have 
an origin in similar unresolved multiple ejections. Detailed 
inspection of radio light curves around these periods may 
shed some light on this. 

The mean position angle of the jet as observed by MER- 
LIN, 142 degrees, is consistent with the range of 140 - 160 
degrees observed with the VLA (RM99). 



4.3 X-ray and radio state at jet formation 

GRS 1915+105 offers a unique opportunity to study the 
jet - disc relation, and in this case we can directly relate 
changes in the accretion disc X-ray emission to the ejection 
of significant quantities of material from the system. 

Fig. 1 shows the evolving sequence of events as observed 
by radio and X-ray monitoring of the source. The onset of 
the plateau phase is preceded by a significant radio flare, 
similar enough to those mapped by us after the plateau state 
to interpret it as a major disc ejection. 

Almost immediately following this ejection the plateau 
state is established. It was already known from PF97 that 
this state, with hard, quasi-steady X-ray emission, corre- 
sponded to bright and fairly steady radio emission at 15 
GHz. Combination of the RT and GBI data now clearly 
show that this state is indeed radio bright, with a flat or in- 
verted spectral index. This is highly indicative of absorbed 
emission (although it is conceivable that the electron ac- 
celeration mechanism changes to produce a different, op- 
tically thin, spectrum which mimics absorption), and may 
suggest the formation of a large optically thick jet. If so, it 
seems possible that this may be associated with the infrared 
jet imaged by Sams, Eckart & Sunyaev (1996). The X-ray 
emission may also support such a picture, as the soft (disc) 
component appears to be very weak, if present at all, dur- 
ing this state (Mendez & Belloni, private communication). 
In addition, no radio oscillations, associated with inner disc 
instabilities, have been observed during the plateau state. 

Fig. 7 shows in more detail the radio and X-ray state 
of the source at the times of formation of the jets we have 
imaged with MERLIN. In particular note the four panels on 
the right hand side of Fig. 7, showing the 15-GHz RT mon- 
itoring of GRS 1915+105 between MJD 50750 and 50754, 
i.e. between the first ejection, corresponding to NW+SE, 
and the second, corresponding to CI. It appears that for 
the entire period between these two ejections, the inner ac- 
cretion disc is unstable and material is repeatedly advected 
and ejected on timescales of tens of minutes (see Belloni 
et al. 1997a,b; Fender et al. 1997; PF97; Eikenberry et al. 
1998; Mirabel et al. 1998; Fender & Pooley 1998). No such 
periods of large amplitude inner disc instabilities are obvi- 
ous between subsequent ejections, although the ejection of 
C2 appears to correspond to a minor X-ray flare, and C3 
possibly corresponds to a brief X-ray dip. 
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Figure 7. A detailed view of the radio and X-ray behaviour of GRS 1915+105 during the ejection episode. Symbols are as shown in 
Fig. 1. The vertical dotted lines, with error bars, indicate our derived times for the beginning of each of the four ejections. The four 
solid vertical bars labelled a — d drawn on the top left radio monitoring panel between the first and second ejections, indicate the four 
periods of 15-GHz monitoring highlighted in the panels on the right. The source appears to exhibit continuous short period (20-40 min) 
oscillations in this phase, indicative of repeated advection or ejection and refill of the inner accretion disc. The radio QPO do not repeat 
after the second ejection, but there are some indicators of the second and third ejections in the X-ray data. 



5 CONCLUSIONS 

Our MERLIN observations have revealed four major rela- 
tivistic ejections from GRS 1915+105 over a period of ap- 
proximately two weeks. Over ten epochs of observation wc 
have consistently measured proper motions of 23.6 + 0.5 and 
10.0 ± 0.5 mas d _1 , for the approaching and receding com- 
ponents respectively. The proper motion of the approach- 
ing component is more than 30% higher than that reported 
from VLA observations of the source (MR94; RM99). Un- 
der the assumption of intrinsic symmetry, we have derived 
an upper limit for the distance to the system of 11.2 ± 0.8 
kpc. While compatible with the distance of 12.5 ± 1.5 kpc 
quoted by MR94, it seems that the favoured distance to 
the system should be revised downwards by at least 1 kpc. 
Solving for angle to the line-of-sight and intrinsic velocity, 
again under the assumption of an intrinsically symmetric 
ejection, we cannot avoid deriving a significantly higher ve- 
locity than MR94, by around A/3 = 0.1. We have investi- 



gated whether or not the lower resolution of the VLA could 
result in lower proper motions being measured, as a result 
of multiple blended components and an increase in decay 
rate with distance from the core, and this does seem possi- 
ble. On the other hand, there is no reason to believe that 
the jet velocity should be fixed - of all the astrophysical ob- 
jects with relativistic jets, only SS 433 (jet velocity 0.26c, 
e.g. Vermeulen et al. 1993) is established to have a constant 
jet velocity. Further observations at high angular resolutions 
with MERLIN and/or VLBI, preferably simultaneous with 
lower resolution VLA observations, will be required to inves- 
tigate whether there is a systematic deceleration of ejecta on 
angular scales of > 0.3 arcsec. If MR94 are correct in infer- 
ring a large distance for GRS 1915+105 from the column 
density to the source, then the bulk velocity of the outflow 
is almost certainly much higher than the currently accepted 
value of 0.92 c; at 11 kpc we derive 0.98t ;^c at 66 ± 2 
degrees to the line-of-sight. 

Our polarisation observations clearly reveal rapid evo- 
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lution of the magnetic field in the ejecta on timescales of a 
day or less. This would seem to imply that the region from 
which the polarised emission arises is smaller than one light 
day across; consideration of time dilation (stretching intrin- 
sic timescales when observed in our frame, given our derived 
solutions for 6 and 0) only makes this size smaller. The de- 
creasing trend of the polarisation with distance from the core 
suggests increasing randomisation of the field as the ejecta 
evolve, although multi-frequency polarisation measurements 
are required to rule out or constrain Faraday rotation effects. 

As already noted by MR94 and others, the power re- 
quired for the formation of the jet is immense, far greater 
than the Eddington luminosity for a solar mass object (at 
11 kpc), even without the inclusion of a proton content. Jet 
formation may well be the dominant power output chan- 
nel during such periods, and possibly also during periods of 
smaller oscillation ejections. In addition, when one proton 
per electron is added the minimum mass flow rate becomes 
comparable to the mass accretion rates derived from X-ray 
spectral fits. So, it is possible that during jet formation pe- 
riods a significant fraction of the inflowing mass is expelled 
and does not fall into the black hole. A similar possibility 
exists for the minor oscillation ejections (Fender & Pooley 
1998). Good wavelength coverage of flares to both higher 
and lower frequencies is required to better determine the lu- 
minosity, energy and mass of the ejections. An accurate si- 
multaneous comparison of mass flow through the disc (from 
X-ray spectral fits) and outflow rates (from radio observa- 
tions) would be of great interest, to see whether the flows 
are advection- or ejection-dominated. 

Radio and X-ray monitoring of GRS 1915+105 (e.g. 
Harmon et al. 1997; PF97) had already hinted at a delayed 
relation between hard X-ray states and radio outbursts. 
These MERLIN observations have established that these ra- 
dio outbursts do indeed correspond to relativistic ejections 
following plateau states. The nature of these states, and the 
associated inverted spectrum radio emission is still unclear, 
and warrants further study. It is also of interest to clarify 
whether or not the radio flare which precedes the plateau 
state (Fig. 1) also, as expected, corresponds to a major ejec- 
tion. Why the accretion disc - jet system in GRS 1915+105 
appears to switch so rapidly between major mass ejections, 
short period (oscillation) instabilities and back again (Fig. 7) 
is still very uncertain. Atoyan & Aharonian (1997) discuss 
the disruption of the accretion disc due to recoil momen- 
tum from an asymmetric ejection. However, it is difficult to 
reconcile their model with the apparent very rapid reforma- 
tion of the (unstable) inner accretion disc after the major 
(NE+SW) ejection. Finally, it is tempting to ascribe the ap- 
parent 'period' of the jet bending of ~ 7 days to an orbital 
period, as was found to be the case for GRO J1655-40, but 
at present this is no more than speculation. 

These MERLIN observations are further evidence that 
GRS 1915+105 repeatedly produces relativistic ejections of 
massive clouds of synchrotron-emitting electrons. Our ob- 
servations of very high proper motions in the inner 0.3 arcsec 
of the jet show that galactic stellar-mass black holes are ca- 
pable of accelerating matter to velocities very close to the 
speed of light. 
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